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معرفی دو اصلاح مو ثر برای بهبود پیش‌بینی رفتار دینامیکی خاکت‌های دانه‌ای در مدل رفتاری 
پاستور - زینکوویچ * 


امین ایرجی( 
مین ابر جی 


چکیده دو اصلاح موثر برای بهبود پیش‌بینی رفتار دینامیکی خاک‌های دانه‌ای توسط مدل رفتاری حمیری تعميم‌یافته )0 
((۵:/>1/ پاستور- زینکوویچ (2ظ ,۵۵0۲-200106۷ ) معرفی شده‌است. اصلاح اول برای در نظ رکرفتن پدیدهُ اثر سست‌شوندگی در 
با رگذاری متناوب عاک های دانه‌ای متراکم و اصلاح دوم برای بهبود پیش‌بینی پد ید رچتینگ (0161618]) یعنی کرنش‌های حمیری 
پیش‌رونده در عاک‌های دانه‌ای تحت بارگذاری متناوب اعمال شده‌است. از سه آزمایش سه‌محوری و سه آزمای شکرنش مسطح تحت 
با رگذاری متناوب برای صحت سنجی مدل اصلاح‌شده استفاده شده‌است. در هر دو اصلاح مذکور. پیش‌بینی تغییرات تنش در مقابل کرنش و 
هم‌چنین تغییرات کرنش حجمی رضایت‌بخش بوده‌اند؛ با این حال نتایج پیش‌بینی تغییرات کرنش حجمی در اصلاح دوم نسبت به اصلاح اول 


تطابق بیشتری با آزمایش‌ها داشتند. 


واژه‌های کلیدی نظریةٌ خمیری تعمیم‌یافته» اصلاح مدل رفتاری. خاک‌های دانه‌ای متراکم» شبیه‌سازی بارگذاری متناوب. سست‌شوندگی: 
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"تاریخ دریافت مقاله ٩۳/۱/۲۲‏ و تاریخ پذیرش آن ۹۷/۱۲/۱۲ می‌باشد. 


(۱) نویسنده مسئول استادیار دانشكده فنی مهندسی خوی, دانشگاه ارومیه. مه هتصت ۵ ززحز یه +اتقصصرز 
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مقد مه 
تحلیل و طراحی دقیق سازه‌های ژئوتکنیکی تحت بارهای 
دینامیکی و لرزه‌ای یکی از چالش‌های مهندسان و محققان 
در دهه‌های اخیر می‌باشد. بسیاری از سازه‌های 
ژئوتکنیکی در شرایط بهره‌برداری عادی و تحت بارهای 
استاتیکی به‌طور معمول معیارهای طراحی را برآورده می- 
کنند. با این حال در شرایط ویژه از جمله بارهای لرزه‌ای 
مشکل پایداری یا تغییرشکل‌های بیش از حد مجاز بروز 
می‌کند که امکان بهره‌برداری را مختل می‌کنند. به‌طور 
کلی روش‌های مذکور را می‌توان در چند دستهُ اصلی 
شامل تعادل حدی. تحلیل حدی و روش‌های عددی 
دسته‌بندی کرد. از آنجا که تعیین مقاومت و تغیبرشکل 
سازه‌ها همواره دو هدف اصلی محفقان در این خصوص 
بوده‌است. لذا روش‌های عددی همواره کامل‌ترین روش 
بهعصوص برای شبیه‌سازی دینامیکی معرفی شده‌است. 
مدل رفتاری دینامیکی مصالح به‌عنوان اصلی‌ترین مسأله 
در روش‌های عددی دینامیکی مطرح می‌باشد؛ خصوصاً 
این که رفتار مصالح ژئوتکنیکی عمدتاً به‌صورت غیرخطی 
و ارتجاعی خمیری می‌باشد. محفقان مختلف مدل‌های 
رفتاری دینامیکی متنوعی را برای پیش‌بینی رفتار مصالح 
ژئوتکنیکی معرفی کرده‌اند. با این حال تعداد بسیاری از 
آنها تنها در حالت‌های مشخص تنش و شرایط مرزی» 
جواب‌های قابل قبولی ارائه می‌دهند؛ لذا اهمیت تحقیق 
و توسعه در خصوص مدل‌های رفتاری دینامیکی بیش از 
پیش نمایان می‌شود. 

مدل رفتاری ۳2 (۳۵۵]۵۳-1601060«1627) بر پاية 
تئوری خمیری تعمیم‌یافته تعریف شده‌است. نظریة اخیر 
توسط مروز و زینکوویچ [1] و زینکوویچ و مروز [2] 
معرفی شده و توسط پژوهشگران دیگر از قبیل زینکوویچ 
و همکاران [3] و پاستور و همکاران [4-6] توسعه داده 
لایر بل تقامی انس با رقار مانت‌های سس 
و متراکم را تحت شرایط استاتیکی و دینامیکی پیش‌بینی 
نماید. در مدل ۳2 از سطوح تسلیم و پتانسیل خمیری 
استفاده نشده و به‌جای آنها از گرادیان‌های سطح تسلیم 


نشريه مهندسی عمران فردوسی 


معرفی دو اصلاح موثر برای بهبود پیش‌بینی رفتار دینامیکی حاک‌های... 


و یتانسیل تممیری اسفاده؛ شده‌است: .این گرادیان‌ها خر 
رابطه بین نرخ تنش و نرخ کرنش لحاظ می‌شوند. 
مهم‌ترین مزیت مدل ۳2 حالت نموی رابطه بین تنش و 
کرنش است که قادر است به‌آسانی رفتار خاک دانه‌ای را 
تحت بارگذاری‌های پیچید؛ دینامیکی شبیه‌سازی کند. 
پاستور و همکاران [6] این مدل را برای شبیه‌سازی 
روان گرایی و رفتار زهکشی‌نشده ماسه به‌کار بردند. 
پژوهشگران مختلفی برروی نظریةٌ خمیری تعمیم- 
یافته مطالعه نموده و آن را به روش‌های مختلفی توسعه 
داده‌اند. . پانتون ی همکاران [7]. اف ناهمسانی 
(۸۹0۵۲۲0۵۲) را در این مدل لحاظ کرده‌اند. باهدا و 
همکاران [8] .عدل اصلی زا تعییر داده‌اند. و اثر پارامتر 
حالت و قانون سخت‌شوندگی دوگانه‌ای ۲۵016) 
(00۵6 عمندعل‌هط را در آن اعمال کرده‌اند. هم‌چنین 
نظریةٌ خمیری تعمیم‌یافته برای مدل‌سازی خاک‌های 
نیمه‌اشباع و با روش عددی انتگراسیون غیرصریح 
(0مصاعجظ صمتاهععط1 ا1۳00110) توسعه داده شده‌است 
[9]. در مطالعة دیگر پدیدهٌ تشکیل باند برشی 50627) 
(0۵00 در خاک‌های کاملاً اشباع و نیمه‌اشباع با روش 
ممتاحعنلمعم1 صتدتاه عتصقطیه و محیط ‏ چندفازی 
(فازهای جامد. مایع و گاز» بااستفاده از نظرية حمیری 
تعمیم‌یافته بررسی شده‌است [10] . لینگ و لیو [11] اثر 
وابستگی پارامترها را به سطح تنش و رفتار تراکمی ماسه 
تحت بارگذاری یک‌سویه و متناوب در نظر گرفتند. 
مرودو و همکاران [12] یک مدل توسعه‌یافتة خمیری 
تعمیم‌یافته را ارائه نمودند که فادر به شبیه‌سازی پدیده 
شکست ۵ در مصالح 
ژئوتکنیکی می‌باشد. میرا و همکاران [13] در بخش 
الاستیک مدل و تغییرشکل‌های برگشت‌پذیر به‌جای 


(مصممصم صمام 


فرمولاسیون هیپو - الاستیک (3(00-012500) که در مدل 
باه کار رده نیاو ترتو تبون هپیرت ایک 
(۲۱۲0۵-۵125010) استفاده کر دند. 

لشکری و لطیفی [14] یک اصلاحیه برای بهبود 
پیش‌بینی پدیدهُ روانگرایی در ماسه توسط مدل رفتاری 
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پاستور- زینکوویچ- چان معرفی کردند. این اصلاحیه 
قابلیت پیش‌بینی روانگرایی ماسه در هنگام چرخحش 
معروهای بای اضلین: را نف تعلالعل کوی اضاقه ی کل 
گلچین و لشکری [15] پدید؛ُ تغییرات ویژگی‌های 
ارتجاعی مصالح دانه‌ای در اثر تغییرات کرنش خمیری را 
با وارد کردن یک مدل هیپرالاستیک در چهارچوب نظرية 
خمیری سطح مرزی حالت حدی بررسی کردند. شبیه- 
سازی‌ها حاکی از رضایت‌بخش بودن نتایج به‌دست‌آمده 
بود. ایرجی و همکاران [16] یک اصلاح موثر برای بهبود 
پیش‌بینی پديدة اثر سست‌شوندگی ماسهٌ متراکم تحت 
بارگذاری متناوب را در چهارچوب نظریة خمیری 
تعمیم‌یافته ارائه کرده‌اند. گورانی و حمیدی [17] یک 
مدل خمیری تعميم‌يافتةُ اصلاح‌شده را برای پیش‌بینی 
رفتار مصالح دانه‌ای حاوی مخلوط ماسه و شن ارائه 
کرده‌اند. آنها برای این منظور خط حالت بحرانی مدل 
مذکور را برای مصالح مخلوط شن و ماسه اصلاح کرده- 
اند. روانبخش و حمیدی [18] با اصلاح خط حالت 
ار 
چهارچوب جدیدی را برای پیش‌بینی رفتار این نوع ماسه 
ارائه کردند. نجمه و لطیفی [19] یک روش تحلیلی برای 
روانگرایی ماس سست و برمبنای مفهوم و تعریف پایه‌ای 
روانگرایی ارائه داده‌اند. آنها برای این منظور از مدل 
رفتاری دافالیاس - منظری (۳1207271-عهن1ه۲(۵1) استفاده 
کرده و نسبت تنش در لحظه شروع روانگرایی را 
به‌صورت تابعی از پارامترهای مدل ارائه داده‌اند. 
ماسه‌های سست و متراکم. طی بارگذاری دینامیکی 
رفتارهای متفاوتی را در باررگذاری مجدد بروز می‌دهند. 
در طول بارگذاری مجدد برروی ماسهٌ سست؛ سختی 
ماس سست افزایش پیدا می‌کند. این پدیده در مدل 
رفتاری ,۳2 لحاظ شده‌است؛ با این حال در ماسة متراکم 
به‌علت دانسیتهٌ نسبی بالا درطی بارگذاری مجدد. کاهش 
سشتی و رفتار اتساعی مشاهده می‌شود؛ که مدل ۳2 قادر 
به شبیه‌سازی چنین رفتاری نیست و یک اصلاح مژثر در 


این خصوص ضروری به‌نظر می‌رسد. از سوی دیگر در 
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مدل ۳2 تغییرشکل‌های نهایی پیش‌بینی‌شده در نمونه- 
های تحت بارگذاری متناوب. به‌ حصوص با سیکل‌های 
کوچک بارگذاری و باربرداری بسیار بیشتر از نتایج واقعی 
می‌باشد. به‌عبارت دیگر مدل ۳2 قادر به پیش‌بینی 
مطلوب پدیدهٌ رچتینگ نیست و یک اصلاح دیگر در اين 
خصوص ضروری به‌نظر می‌رسد (لازم به ذکر است که 
مفهوم پدیدهٌ رچتینگ بیشتر در مهندسی مکانیک و 
متالورژی کاربرد دارد و در اين رشته‌ها نیز تاکنون معادل 
خاس «فاشیم. بای ر اف کلجه ات تفنواست؛ 
بیشترین کاربرد این واژه در ادبیات فنی به‌صورت 
«رچتینگ) است. با این حال در موارد کمی از عبارت 
«حزش سیکلی) نیز استفاده شده‌است). اصلاحيهٌ اول 
مدل منجر به افزایش تعداد پارامترهای مدل نمی‌شود 
درحالی‌که در اصلاحيهٌ دوم تعداد پارامترهای مدل 
ه‌اندازٌ دو پارامتر افزايش می‌یابد. سه آزمایش سه- 
محوری و سه آزمایش کرنش مسطح تحت بارگذاری 
متناوب برروی خاک‌های ماسه‌ای برای صحت‌سنجی 
استفاده شده‌اند. تغییرات تنش درمقابل کرنش و هم‌چنین 
تشییرانت: کرفیی اجبی تبزرسین قبله اب و کمام شییات 
سازی‌های انجام‌شده شرایط زهکشی‌شده در نظر گرفته 
شده‌است. خاطرنشان می‌شود که مدل اصلاح‌شده تنها 
توانایی شبیه‌سازی رفتار دینامیکی در شرایط زهکشی‌شده 
را دارد و قابلیت شبیه‌سازی در شرایط زهکشی‌نشده را 
ادن خویای داتسرف 
((۳0۳۸) و سپس بااستفاده از زبان برنامه‌نویسی 


فیش (۳15۴8) در نرم‌افزار ۳1,۸6 نوشته شده‌است. 


تشریح مدل پاستور- زینکوویچ 
مدل ۱۳2 ساسا در ف ضای سه‌محوری 8,0,6 تعریف 
شده‌است. این مدل در بارگذاری یکسویه (۷]۵۵06010) 
به ۷ پارامتر و در بار گذاری متناوب به ۱۰ پارامتر نیاز 
دارد. 
رابطة بین نرخ تنش و نرخ کرنش برای مص‌الحی 
که رفتار آنها تابع زمان نیست. به‌صورت زیر تعریف 
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می‌شود: 

)۱ 56 < 06 
که در آن 0 و ۲۳ به‌ترتیب نرخ تانسور تتش و 

تانسور سختی الاستوپلاستیک هستند. عملگرد «» بیانگر 

ضرب عددی دو تانسور است. تانسور سختی 

لاستربلازشیگ »در رنه ری تم باه با مروت 

زیر تعریف می‌شود: 


هنن ررمهنگظ 


0۶ - ۶ (۳ 


هن ۲ 

که در آن ۰۳۴ ظ روط ون/ 11 به ترتیب تاذ سور 
سختی الاسستیک. بردار جهت بارگذاری و باربرداری» 
بردار جهت جریان خمیری تحت شرایط بارگذاری یا 
باربرداری و ضریب پلاستیک بارگذاری يا باربرداری 
می‌باشند. هم‌چنین اندیس‌های بآ و ۲ به‌ترتیب حالت 
بارگذاری و باربرداری را نشان می‌دهند. 

رفتار الاستیک از قانون هوک (1060) پیروی 
می‌کند. ضریب برشی ارتجاعی توسط رابطهٌ زیر تعریف 
می‌شود: 
۳ (ش) ,6 < 6 

در رابط؛ بالا 9 تنش مژثر میانگین 
3 < ۰۳ و۲ فشار اتمسفر, رآ اولین نامتغیر تتش و و 
مدول پرشی پایه می باش‌ند. هم‌چنین مقدار 0 
۶۵ تحاظ می‌شود. مدول حجمی ارتجاعی 


به‌صورت زیر تعریف می‌شود: 


1 - )( 3 


3 
که در آن 16 مدول حجمی پایه است. 
براساس رابطة پیشنهادی نووا و وود [120 
نمودند که در مدل ۳2 به‌کار برده شده‌است: 
۹ 
- 2۷ 
(۵( (- 04( + ) << یا 
که در آن,غ و جع به‌ترئیب نرخ کرنش حجمی 


خمیری و نرخ کرنش انحرافی خمیری: با شسیب حط 
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معرفی دو اصلاح موثر برای بهبود پیش‌بینی رفتار دینامیکی حاک‌های... 


حالت بحرانی در ف ضای 0 - ۲ 0/۲ < 1 ذسبت تنش 
0 با رامع خانت سل هس نیلوا واوية ام طکاک 
داخلی در حالت بحرانیع و زاو بة لود (1,000) 0 


روابط زیر را دارد: 


65۹116 
)3 661030 51۳00- 3 ٍِ ۷ 
"( قط تد < 30وزو 

2 3 

12 

که در آن و[ و و[ به‌ترتیب نامتغیرهای دوم و سوم 
بردار گرادیان سطح 7 سلیم بردار جهت بارگذاری 
یعنی و بردار گرادیان سطح پتانسیل خمیری بردار 


خمیری (۲۷16 ۳۱0۷ غیروابسته است. به‌طوری که دو 
بردار گرادیان ذکر شده یک‌سان نی ستند و به صورت زیر 
تعریف می‌شوند؛ 


0۳0۵ ۹99 


۳ 0 


< )06 1, - 002605 30 /2(* 
(۸) 


0 08 "ِ 
۳ 


0۳999 


2(۳/ 0530 و/(9- ,1 ,ع0) < 


(8) 


]و 8 به‌ترتیب سطوح ‏ سلیم و پتاا سیل خمیری 
مدل هستند. 1/۶ ثابت مدل است که در اندازه سطح 


تسلیم تأثیرگذار است و (1- 00۶ (0 + 1) < مه. 
ضریب خمیری در شرایط بار گذاری اوایه 

به‌صورت زیر تعریف می‌شود: 

(0۰ (,۱ + ,)ما > را 
0 ضریب خمیری اولیه و م11 ,لآ و ول ضرایب 

خمیری هستند که به‌صورت زیر تعریف می‌شوند: 


۳ 
11 - 1 ِ ۱۱) 
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امین ایرجی 
که در آن: 
۱ ,۷( + 1) < 7 
۱۳ 7 1 2 پآ[ 
۱ (۵8-) 6086۵0 < وا 


۶ پارامتر نسبت تنش است.80 و 8 پارامترهای 
ثابت مدل ه ستند. همچنین 5 کرنش انحرافی پلا ستیک 
تجمعی اسست: |غ| ل < و غ نرخ کرنش انحرافی 
خمیری است. 

ضریب خمیری برای بارگذاری مجدد به‌صورت 
زیر تعریف می‌شود: 

(۱۵( پردلا(ی + واامتامل! ع رل 


در رابطة بالا ضریب »,رو به‌منظور در نظر گرفتن 


تاریخچه تنش می‌باشد و به‌صورت زیر تعریف می‌شود: 


6۳ "سس ررو11 
۱۳( (گ -1). < 


در رابطة بالا ) تنش بسیج‌شده 4عتنان0/0) 
(۳655ا٩‏ در هر لحظهٌ باررگذاری آانبتتتا: هم‌چنین / پارامتر 
ثابت مدل رفتاری می‌باشد. ضریب خمیری مدل در 


باربرداری به‌صورت زیر تعریف شده‌است: 


۱ ۷ 
101 )| م۲ كت ۲ 
)1 


من < ول 


)۱۸( 
در اين رابطه. ون[ ضریب خمیری باربرداری اولیه 
می‌باشد. ,[1 نیز نسبت تنش در نقطه‌ای است که از آنجا 
باربرداری شروع می‌شسود. ن۷ یک پارامتر ثابت مدل 

است. اندیس آ] حالت باربرداری را نشان می‌دهد. 
بردار پم در هنگام باربرداری به‌صورت زیر تعریف 
شده‌است: 


(رروت رود ,|رجت|-) ۳ (ویظ ,رو ریص) 2 بع(1 
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۷۷ 


(۱۹( 
بردار 8 در بارگذاری و باربرداری یکسان است؛ 


بنابراین داریم: 


۲( (یدره > (نه ندرم > ده 


تعیین پارامترهای مدل ,۳7 
در مدل رفتاری ۰۳2 مت از شیب اولیةٌ نمودار تنش- 
گرنکن تهاستا, ی آبله و1 ان آزهایکن فشاری. مک 
دما 60۳50۲6951010 150170016) يا با برازش دادن شیب- 
های اولی نمودار ل - رع یا و - رع در آزمایش و شبیه- 
سازی حاصل می‌شود. م1 شیب نمودار 8 - 0 در حالت 
بحرانی است. برای تعیین مقدار اولیه برای 1۶ معمولاً از 
فرمول تقریبی «آحی۱۷]:۷ استفاده می‌شود. این نسبت 
برای ماسة متراکم حدود ۰/٩‏ الی ۱/۱ و برای ماسة سست 
حدود ۰/۳ الی ۰/۳ است. مقدار تقریبی 16۶ از ضرب این 
نسبت در مقدار م1 به‌دست می‌آید. مقدار دقیق آن نیز از 
کالیبره کردن آن با نتایج آزمایشگاهی تعیین می‌شود. 
مقدار » از شیب نمودار اتساع () - نسبت تنش (0) 
به‌دست می‌آید (رابطٌ ۵). با این حال این پارامتر در بیشتر 
شبیه‌سازی‌ها ۰/4۵ فرض می‌شود. ول[ در آزمایش سه- 
محوری يا کرنش مسطح از شیب اوليةٌ نمودارهای تنش 
پرشی- کرنش محوری يا شیب اولیةٌ نمودار کرنش 
حجمی- کرنش محوری در آزمایش‌های زهکشی‌شده یا 
شیب اولیةٌ مسیر تنش "0 - 0 در آزمایش‌های زهکشی- 
نشده به‌دست می‌آید. 80 و ,8 از نمودارهای تنش برشی- 
کرنش محوری آزمایش به‌دست می‌آید. مقادیر 80 و 8۰ 
بااستفاده از روش سعی و خطا و طوری انتخاب می‌شوند 
تا نمودار شبیه‌سازی شده بیشترین تطابق را با نمودار 
آزمایش داشته باشد. مقادیر اولیةٌ توصیه‌شده برای این 
پارامترها به ترتیب ۶/۲ و ۰/۲ است. بازُ توصیه‌شده برای 
این پارامترها به‌ترتیب برابر با ۱/۵ الی ۵/۰ و ۰/۱ الی ۰/۲ 


هستد. 
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۷۸ 


من[ از شیب اولية نمودار تتش برشی - کرنش 
محوری در اولین باربرداری به‌دست می‌آید. ,۷ از 
کالیبره کردن با نرخ تغییر شیب نمودار تتش برشسی- 
کرنش محوری در اولین باربرداری به‌د ست می‌آید. ۷ از 
کالیبره کردن با شیب نمودار تتش برشسی- کرنش 
محوری در اولین بارگذاری مجدد به‌دست می‌آید. واحد 
0 و م16 برحسب کیلوپاس کال و بقيةٌ پارامترها بدون 


اصلاح پیشنهادی اول: اعمال اثر سست‌شوندگی 
پارامتر مرول در مدل ۳2 اثر تاریخ چذه تنش را در 
بارگذاری مجدد و در بارگذاری متناوب نشان می‌دهد 
۱ ترا این تاش سم صرق از گر 
به‌عبارت دیگر کاهش سختی در بارگذاری مجدد ما سة 
متراکم 1حاظ شود م قدار مرو[ بایستی در بازه 
| 
ابیت که قدای ون هقی نا شاه تا که 
اه مان کر مان رم ارم هت ای یط 
مقدار بیشینه کمیت ؟ در پایان بارگذاری اولیه ذخیره 
می‌شود و برای محاسبهة مرول در بارگذاری مجدد 
استفاده می‌شود. مقدار مرول در بارگذاری اولیه برابر 
یک است. در ابتدای بارگذاری همواره شرایط برومیة > ؟ 
برقرار است. بنابراین مطابق رابطة (۱1) مقدار مروا از 
یک بیشتر می‌شود. با این حال زمانی که مقدار 6 در طول 
بارگذاری افزایش می‌یابد و از رو نیز بٍ شتر می شود 
ی 
مقدار برواا دربرابر کرنش برشی در اولین بار گذاری 
مجدد در یکی از آز مایش های سه‌محوری دینامیکی 
زهکشی‌شدة انجام شده توسط پرادهان و همکاران [21] 
برروی ما سه متراکم و در حالت شش (200605100) را 


نشريه مهندسی عمران فردوسی 


معرفی دو اصلاح موثر برای بهبود پیش‌بینی رفتار دینامیکی حاک‌های... 


۱1 


04 0.02 0 0002 
و۱ 


شکل ۱ تغییرات مرو در برابر کرنش برشی در مرحلهٌ بارگذاری 


مجدد یک آزمایش سه‌محوری متناوب (پرادهان و همکاران [21]) 


برروی ماسةٌ متراکم در مدل ۳72 و در حالت کششی (0625100ه) 


تخاک‌های شاجم‌ای مشسست یا تضتعا . سب ده 
بارگذاری‌های مجدد در طول بارگذاری متناوب. رفتاری 
سخت‌شونده دارند. بنابراین می‌توان نتیجه گرفت که 
نحوة تعریف پارامتر م11 در مدل ۳2 مناسب خاک‌های 
دانه‌ای سست است. در ماسه‌های متراکم همان‌گونه که 
نتایج آزمایش‌های پرادهان و همکاران [21] و ماسودا و 
همکاران [22] نشان می‌دهند. نمونه در بارگذاری مجدد 
کاهش سختی دارده شیب منحنی بارگذاری آن کمتر است 
و تراکم نمونه‌ها کمتر می‌شود. بنابراین پارامتر مرو در 
مدل ۳2 قادر به شبیه‌سازی چنین رفتاری نخواهد بود. 

پرای شبیه‌سازی نرم‌شوندگی ماسة متراکم لازم 
است تا ضریب ,رو مقداری بین صفر و یک داشته باشد. 
هم‌چنین مقدار ۷ باید منفی باشد. بنابراین در رابطة (۱1) 
بایستی یک عبارت در 7 ضرب شود تا علامت آن منفی 
شود. این عبارت بایستی طوری تعریف شود که قادر به 
تشخیص حالت متراکم يا سست خاک دانه‌ای باشد. در 
تعریف جدید به‌جای ۷ از 7 استفاده می‌کنيم. رابطه 
جدید زیر به‌عنوان جایگزین رابطةٌ (۱3) معرفی 
شده‌است: 


وم 
و7 05 ۳ 


(۲۱( تس < مرول 


5 
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امین ایرجی 


غبارت 2 5 حالت سست يا متراکم خاک 
دانه‌ای را تشخیص می‌دهد. بازهٌ معمول برای دانسیتة 
نسبی خاک‌های داله‌ای بیرن ۰/۳ 9 ۰/۷ است. بثایراین 
مقدار ۰/۵ به عنوان دانسیتة نسبی متوسط در نظر گرفته 
می‌شود. مقدار ی نیز تقریباً مساوی با دانسيتة نسبی 
خاک () در نظر گرفته می‌شود. مطابق رابطه (۲۱) 
هنگامی‌که نوع خاک ماسة متراکم است. عبارت 
2 0( قی ی و ترس وه این که مار 
شتا اهزاین فان زوا ها از بخ‌سرا هد 
بود. با ضرب مقدار »,و11 در ,رل ضریب خمیری در 
بارگذاری مجدد کاهش خواهد یافت. هنگامی که نوع 
فاگ فا نگ ارت ادا ( - 05) ۱۹ 
بنابراین مقدار مرو برابر يا بزرگ‌تر از یک خواهد شد 
که نتیجة آن بزرگ‌تر بودن سختی ماسةٌ سست در 
با گذار ی تمه اعد شه, معادی اباراش‌ها وا می نان 
طوری کالیبره کرد که مقدار عبارت هم - 605 در 
مدل ۳2 اصلاح‌شده با پارامتر 7 در مدل ۳2 در ماسة 
سست یکسان شود. بنابراین اصلاح معرفی‌شده تنها در 
خصوص پیش‌بینی رفتار دینامیکی ماسة متراکم کارساز 
است. 

خاطرنشان می‌شود که در صورت تخصیص مقدار 
8 تِ تا 7 ک یا اناوت 
صفر می‌شود. درنتیجه مرو[ در رابطٌ (۲۱) برابر با یک 
خواهد شد. به عبارت دیگر نمونه تراکم متوسطی خواهد 
دز ات رال واط ی یرای ی سوه 
بارگذاری‌های مجدد از بین خواهد رفت. 

در مدل ۳2 اصلاح شده در ابتدای بارگذاری مجدد 
ماس متراکم» مقدار تنش بسیج‌شده () کوچک‌تر از 
مقدار بيشينةٌ آن است. بنابراین مقدار مرول عددی بین 
صفر و یک خواهد بود. زمانی که مقدار تنش بسیج‌شده به 
مقدار بيشينةٌ قبلی می‌رسد و از آن مقدار نیز تجاوز می- 
کند. مروآ! نیز بیشتر می‌شود و به مقدار بیشینه‌اش» یعنی 
یک می‌رسد. شکل (۲) تغییرات مقدار مرآ اصلاح‌شده 
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۷۹ 


را در برابر کرنش برشی در اولین بارگذاری مجدد در 
یکی از آزمایش‌های سه‌محوری دینامیکی زهکشی‌شده 
انجام‌شده توسط پرادهان و همکاران [21] برروی ماس 


متراکم و در حالت کشش (616051070) نشان می‌دهد. 


۱112 


0002 0 0.02 0,04 


صاصاو تحعط۹ 


شکل ۲ تغییرات برو3[ اصلاح‌شده دربرابر کرنش برشی در مرحلهً 
بارگذاری مجدد یک آزمایش سه‌محوری دینامیکی (پرادهان و 
همکاران [21]) برروی ماسة متراکم در مدل ,7و در حالت کشش 


(طمتعصعه) 


اعتبارسنجی اصلاح پيشنهادشده برای اعمال اثر 
سست شوند گی 
دو آزمایش سه‌محوری دینامیکی زهکشی‌شده توسط 
تاتسوکا و ایشیهارا [23) برروی ماس سست و توسط 
پرادهان و همکاران [21] برروی ماس متراکم و دو 
آزمایش کرنش مسطح دینامیکی زهکشی‌شده توسط 
ماسودا و همکاران [22] برروی ماس متراکم توسط مدل- 
های ۳۶ و 12 اصللاخح‌شده شبیه‌سازی شده‌اند. شبیه‌سازی 
ماسة سست توسط مدل‌های ۲2 و ۳2 اصلاح‌شده نتایج 
یکسانی را می‌دهند. با این حال» مدل ۲2 اصلاح‌شده 
نتایج بسیار بهتری را برای ماس متراکم نسبت به مدل ۳2 


ارائه می‌دهد. 
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ماس سست 
شبیه‌سازی آزمایش تاتسوکا و ایشیهارا 

تاتسوکا و ایشیهارا [23] یک سری آزمایش سه- 
محوری دینامیکی برروی ماسة فوجی ریور 2۷۲6 نزدا۳) 
با دانسیته‌های نسبی مختلف انجام داده‌اند. تخلخل این 
نوع ماسه در بازٌ ۰/۵۳ تا ۱/۰۸ متغیر بوده‌است. نمونه‌ها 
اشباع‌شده بوده و ارتفاع انها جت ۱۰و قطر انها جه ۵ 
بوده‌است. 

آزمایش دینامیکی برروی نمونةٌ نسبتاً سست ماسه 
با تخلخل اولية ۰/۷۶ و تنش جانبی 1۳۵ ۲۰۰ انجام 
شده‌است. تنش‌ها و کرنش‌ها بدین شرح تعریف شده‌اند: 
2 تنش محوری مژثر 0 تنش شعاعی مژثر.- 6 <)9 
(6 تنش انحرافی؛ [207(/3 + 6) <)0 تنش میانگین» 
کرنش محوری .۶ کرنش شعاعی» و (۶ - بة ۷62 
کرنش برشی. 

شکل‌های (۳-لف) و (-لف) به‌ترتیب نمودار 
نسبت تنش- کرنش برشی و نسبت تنش- کرنش حجمی 
را در آزمایش نشان می‌دهند و شکل‌های (۳-سب) و (4- 
ب) شبیه‌سازی آنها را بدون اعمال اصلاح مدل نشان می- 
دهند. بررسی نمودارها مشخص می‌کند که مدل ,۲2 رفتار 
دینامیکی ماسه را به‌خوبی پیش‌بینی کرده‌است. همان‌گونه 
که در شکل (۲سب) مشاهده می‌شود. شیب‌های 
بارگذاری, باربرداری و بارگذاری مجدد در شبیه‌سازی 


تطابق خوبی با شیب‌های متناظر در آزمایش دارند. 


له ۵۵6 ۲۳1۵1۵ 
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نها محعم5 


(الف) 
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معرفی دو اصلاح موثر برای بهبود پیش‌بینی رفتار دینامیکی حاک‌های... 


آهمنهلبهنگ 


003 002 001 0 1,- 02.- 
صنحاه تحعط۹ 


(ب) 
شکل ۳ مقايسهة نتایج شبیه‌سازی با نتایج آزمایش سه‌محوری 
زهکشی‌شده دینامیکی ماسةٌ سست. نمودار نسبت تنش- کرنش 
برشی (الف) آزمایش تاتسوکا و ایشیهارا [23] (ب) شبیه‌سازی با 
مدل ۳2 و ۳2 اصلاح‌شده (هر دو یکسان می‌باشند) 


09045 


(۱974 ,۵1 مامتواع۲) صعهعی۲۳ 


۱ 


1 05 0 05 1 
90 


(الف) 


‌تامایسن؟: 


۱ 


(ب) 


شکل ۶ مقايسهة نتایج شبیه‌سازی با نتایج آزمایش سه‌محوری 
زهکشی‌شده دینامیکی ماسة سست. نمودار نسبت تنش- کرنش 
حجمی (الف) آزمایش تاتسوکا و ایشیهارا [23] (ب) شبیه‌سازی با 


مدل ۲2 
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پات سار نی فا سیر رای ماه 
۶/۷ برابر با ۰/۵۳ می‌باشد که از دیدگاه اصلاح 
معرفی‌شده حاکی از یک نمونة تسبتا سست است. مقدار 
توان ,۲10 یعنی ۷ در مدل ۳2 (به‌طوری که نمودار پیش - 
بینی» بهترین مطابقت را با آزمایش داشته باشد) برابر با 
یک است. این توان در مدل اصلاح‌شده برابر با 
۷ - (2 - 05) وروی با افوضین فقلان 1۱6 
تراغ دار 10 0 دا اهاز 
سیکل بارگذاری مجدد بهترین مطابقت را با آزمایش 
خواهد داشت. بنابراین اصلاح ارائه‌شده برای نمونه‌های 
با مقدار ع/:۷ کوچک‌تر از ۰/۵ يا برای ماسه‌های 
سست کاربردی نخواهد داشت. پارامترهای مدل ۳۸2 برای 
ماسهٌ فوجی ریور (1:762 [[د؟) در جدول (۱) آمده‌است 
که از کالیبراسیون با نتایج آزمایش تاتسوکا و ایشیهارا 


[ 22 به‌کست. آملواسستا 


جدول ۱ پارامترهای مدل ۳2 و در ازمایش سه‌محوری دینامیکی 


(تاتسوکا و ایشیهارا [23]) برای ماسهٌ فوجی ریور ۲۷6 زدا۳) 


پارامتر مدل مقدار 
مدول حجمی پایه (6۳2) ,1 ۱۳۹ 
مدول برشی پایه («0۳ 00 :1-۹ 
پارامتر ثابت مدل :11 ۰/۹ 
شیب خط حالت بحرانی پل ۱/۶ 
ضریب خمیری اولیه (بارگذاری) م11 و 
پارامتر ثابت مدل (باربرداری) ون[ 
پارامتر ثابت مدل (بارگذاری مجدد) ۷ ۱/۰ 
پارامتر ثابت مدل (باربرداری) ن۷ ۱/۰ 
پارامتر ثابت مدل 80 5/۰ 
پارامتر ثابت مدل .8 ۰/۲ 


ماسة متراکم 
شبیه‌سازی آزمایش پرادهان و همکاران 
پرادهان و همکاران [21] یک سری آزمایش‌های سه- 


محوری دینامیکی را روی نمونه‌های ماس اشباع تویورا 
(10(007۵) تحت شرایط زهکشی‌شده انجام داده‌اند. 
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۸ 


نمونه‌ها در ابتدای آزمایش‌ها تنها در یک امتداد تحکیم 

شده‌اند. ماسة تویورا متشکل از دانه‌های تيزگوشه تا نسبتاً 

تيزگوشة کوارتز )٩۰/(‏ و دانه‌های سیلیسی (/4) است. 
پارامترهای فیزیکی این ماسه به‌شرح زیر است: 

,146 < لا 0.169 < ووط , 2.64 < 6 

5 2 ونم ,0.977 < بجسه 


تنش همه‌جانبة (0) ثابت در ازمایش‌های سه- 
محوری انجام‌شده 6۳۵ ۹۸/۱ می‌باشد. تنش‌ها و کرنش‌ها 
بدین شرح تعریف شده‌اند: 02 تنش محوری مژثر 67 
تنش جانبی مژث. (ع6-:6 462 تنش انحرافی, 
(20:(/3 + 05) )0 تنش میانگین مزن ,ع کرنش 
محوری. ,۶ کرنش جانبی. و ( - بع )۷ کرنش برشی 
می‌باشد. در آزمایش انتخاب‌شده برای مطالعهً حاضر 
(نمونة شمارءٌ 0۷0۳09 پرادهان و همکاران [21]) 
تخلخل ماسه در تنش همه‌جانبه ۰/۳۵/۳۴ برابر < ه 
3 است که حاکی از تراکم بالای ماسه می‌باشد. 
پارامترهای این ماسه برای مدل ۳2 و ۳2 اصلاح‌شده در 
جدول (۲) معرفی شده‌اند. که از کالیبراسیون با نتایج 
آزمایش پرادهان و همکاران [21] به‌دست آمده‌اند. 
جدول ۲ پارامترهای مدل ۳2 و ۳2 اصلاح‌شده در آزمایش سه- 


محوری دینامیکی (پرادهان و همکاران [21]) برای ماسه تویورا 


پارامتر مدل مقدار 
مدول حجمی پایه (۳2) مک ۳ 
مدول برشی پایه (6۳2) 00 و 
پارامتر ثابت مدل 11۶ ۰/۸ 
شیب خط حالت بحرانی »14 ۱/۵ 
ضریب خمیری اولیه (بارگذاری) و1 ۱۳/۰۹ 
پارامتر ثابت مدل (باربرداری) 110 0 
پارامتر ثابت مدل (بارگذاری مجدد) *۷ ۱/۰ 
پارامتر ثابت مدل (بارگذاری مجدد) ۷۲ ۳/۹ 
پارامتر ثابت مدل (باربرداری)ن ۷ 1۳۳ 
پارامتر ثابت مدل 80 ۷/۰ 
پارامتر ثابت مدل ,8 ۰/۵ 


* فقط فلل 7 فقط مدل ۳7 اصلاح شده 


نشریه مهندسی عمران فردوسی 


۸۲ 


شکل‌های (۵) و (1) آزمایش دینامیکی ماسهٌ متراکم 
تویورا و شبیه‌سازی آن را نشان می‌دهند. مقدار »1۷1۷1 
بیشتر از ۰/۵ یعنی ۰/۵۶ می‌باشد که نشان‌دهنده متراکم 
بودن ماسه امستتا: همان‌طور که مشاهده می‌شود با اصلاح 
انجام‌شده نمودار شبیه‌سازی در شکل‌های (۵) و (1) به 


نتایج آزمایشگاهی نزدیک‌تر شده‌است. 


ادها عناهن افمته1 


04 -0.2 0 002 004 
5۳09۲ ۵ 


(الف) 


۱۱۵۵۵ 2( ط ممناعلنم5 
اعهصه 2 نزمه جوا رازه 


+004 002 0 002 124 
تاد عححط 


(ب) 
شکل ۵ مقايسة نتایج شبیه‌سازی با مدل ۰۳2 ,۳2 اصلاح شده و 
نتایج آزمایشگاهی برای ماسهٌ متراکم (الف) نسبت تنش- کرنش 


برشیء آزمایش پرادهان و همکاران [21] (ب) نسبت تنش- کرنش 


برشی. شبیه‌سازی 


نشریه مهندسی عمران فردوسی 


معرفی دو اصلاح موثر برای بهبود پیش‌بینی رفتار دینامیکی حاک‌های... 


(۱۵۵۵ .اه او مماه اصعصزه ما 


صجتاه ماه ۷ 


(الف) 


۱۱۵0۵7 ۳7 وه صهاتدآصات۹ 


۱ 


1۲ 


(ب) 


۵۵۱ ۳۶ 0قن0هن ۲ ممتاداصنق 


۱ 


شکل 1 مقایسة نتایج شبیه‌سازی با مل ,۳ ۳7 اصلاح‌شده و 


نتایج آزمایشگاهی برای ماسهٌ متراکم» (الف) نسبت تنش- کرنش 


حجمی: آزمایش پرادهان و همکاران [21] (ب) نسبت تنش- 


کرنش حجمی, شبیه‌سازی 
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امین ایرجی 


شبیه‌سازی آزمایش‌های ماسودا و همکاران 
ماسودا و همکاران [22] رفتار تناوبی ماسهٌ تویورا را در 
حالت متراکم در دستگاه آزمایش کرنش مسطح مورد 
آزمایش قرار داده‌اند. پارامترهای فیزیکی این ماسه به‌شرح 


6 2 لا بط 0.162 < مو9 , 2.64 -< .6 
12 2 «زم6 ,0973 < بجره 


آنها تعدادی آزمایش دینامیکی با مقدار فشار جانبی 
ثابت 6۷ انجام داده‌اند. تعدادی از نمونه‌ها در ابتدای 
آزمایش به‌صورت همگن و تحت شرایط < ,6 < »6 
52 78.5 و تعدادی دیگر به‌صورت غیرهمگن و تحت 
شرایط 78.512 < م6 و 29.512 < ,6 تحکیم 
یافته‌اند. ,ب6 تنش عمودی وارد بر نمونه است. تخلخل 
اولیةٌ نمونة تحکیم‌یافته ب‌صورت همگن 0.654 < و6 و 
تخلخل اولیهٌ نمونهٌ تحکیم‌یافته به‌صورت غیرهمگن 
9 < و می‌باشد. ارتفاع نمونه‌ها ه ۲۰ طول آنها 


۱۲ و عرض آنها (در امتداد ,و6) عه ۸ است. نتایج 


حدول ۳ پارامترهای مدل 7 و ۳7 اصلاح‌شده در 


پارامتر مدل 


مدول حجمی پایه (۳2) م16 
مدول برشی پایه (۴0) ون) 
پارامتر ثابت مدل 1۶ 
شیب خط حالت بحرانی و 
ضریب خمیری اولیه (بارگذاری) 11 
پارامتر ثابت مدل (باربرداری) ون[ 


پارامتر ثابت مدل (بارگذاری مجدد) *۷ 


پارامتر ثابت مدل (بارگذاری مجدد) ۷۲ 


پارامتر ثابت مدل (باربرداری)ی ۷ 
پارامتر ثابت مدل 80 


بازامع یت سدل0 


* فقط مدل ۳2 "فقط مدل ۳2 اصلاح‌شده 
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۸۳ 


آزمایش‌ها برحسب زاويةٌ اصطکاک داخلی بسیج‌شده 


67-0 


< و(ومم ۳ و کرنش برشی و - رع ارائه 


00 
شده‌اند. 


در تحقیق حاضر ازمایش‌های شمارءٌ 5-6 6286 و 
2-5 6256 انتخاب شده و شبیه‌سازی شده‌اند. در طول 
مدت بارگذاری تنش جانبی ثابت 1۳2 ۷۸/۵ برای حالت 
تحکیم يافتة همگن (5-6 6856) و تنش جانبی ثابت 1۳2 
۵ برای حالت تحکیم‌یافتةُ غیرهمگن (2-5 6۵96) بر 
مرزهای جانبی مدل اعمال شده‌است. جدول (۳) 
پارامترهای مدل ۳2 و ۳2 اصلاح‌شده را برای ماسة 
یراگن آنماتتی گر تن مسطح نشان می‌دهد. 

پارامترهای مدل ۳2 اصلاح‌شده همانند پارامترهای 
بل هر اون (۳) نارای ناش ۶ 
استفاده خواهد شد. شکل‌های (۸ تا (۰) نتایج 


آزمایشگاهی نمونه‌های 5-6 6886 و 2-5 6286 به‌همراه 


نتایج پیش بینی مدل ۳2 و مدل ۳2 اصلاح‌شده وا تیان 


می‌دهند. 
آزمایش کرنش مسطح دینامیکی برای ماسة تویورا 
مقدار(5-6 0286) | مقدار (2-5 0۵86) 
12-۱۱۰۶ :2-۱-۰ 
کر 2 -2(۳(3۰ 
۰/۹ ۰/۹۵ 
۱/۲ ۱/۱ 
2 2 
۱/۱ ۱۳/۹۹ 
۱/۰ ۱/۰ 
۳/۰ ۱/۰ 
۱/۰ ۱9/۰ 
1/۰ ۱۰/۰ 
۰/۸ ۱/۱ 


نشریه مهندسی عمران فردوسی 


۸ 


شکل‌های (۷-الف) و (٩-الف)‏ ذشان می‌دهند که 
در طول آزمایش دینامیکی سختی نمونه به‌تدریج کاهش 
می یابد. پیش‌بینی مدل‌های آزمایش توسط مدل ۳2 
اصلاح شده نیز ذشان‌دهنده کاهش سختی در سیکل‌های 
بارگذاری مجدد است. درحالی که این کاهش سختی در 
نتایج مدل ,۳2 قابل مشاهده نیست. 

شکل‌های (۸) و (۱۰) نمودارهای کرنش جانبی- 
کرنش قائم را در آز مایش ها و مدل های ,۳2 و ۳2 
اصلاح شده نشان می‌دهند. نتایج حاکی از بهبود نتایج در 
اثر اصلاح اعمال‌شده هستند. 


اج اجه ۲۲1۵۷۲۸ 


(ح .و۲۲6 مصعنو )نز ط.مصعز ۱-6 مصعنق)- را )من 


0005 0 
حادء تفعاه 


5 


(الف) 


ما۱۵ ۶ ات شمه بط صفتااننه 


(طمسوزه ۱ ۲.مصعن۹ )۱( مصع-,مصونو) رنطسنه 


0.05 


0.005 ۵ 
اما کفت‌ ٩‏ 


(ب) 


نشریه مهندسی عمران فردوسی 


ال۱۵ ۱۶ بط حصتافا تاه 


(ط,مسصونه ۱ ۲.محصونه )۸( مصعز٩-مصونه)‏ رتطگ‌سنو 


۵ 0005 


تاد تف‌ط 


-5 


ج 


شکل ۷ مقایسه نتایج آزمایشگاهی کرنش مسطح دینامیکی 
زهکشی‌شدة ماس متراکم با نتایج مدل ۳2 و مدل ۳2 اصلاح شده 
5-6 ۵96 تنش بسیج‌شده- کرنش برشی الف) آزمایش ماسودا و 
همکاران [22]. ب) شبیه‌سازی با مدل 2 اصلاح‌شده و ج) شبیه - 


سازی با مدل ,۳7 


0 ۲۵۷/6۷۴۵" 
سا :هاعاد ادنانما 
752 
7852 ۳۹ 
مه 0125 < رو ك 
۳ 
ادها هار۵ ۲-1 
(5-6 ودع6) ۵ 
موه 5 
5 
ت 
[3) 5 ۱۷۵۲۷۵۵1 
(الف) 
6 
0.5 
04 
03 ند 
رم 2 
۳ 
5 
1 - 
كِ 
5 
ت‌ِ ‏ 
41 3 
032 
03 


8 ا ا ‏ ها ا ‏ اس( 
(9۵) واه اهنا۱۷۵۲ 


(ب) 
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< 8 08 

۵۵1 ۳7 بح ( ما ات٩‏ سس ِ 

3 05 

۳ 

3 04 

ی 

سته رم 

-- ِِ 

2 8 02 

و ِ 

۳ 4 

ِ 8 ۰ 

۳ ۹ 

02 

ِ 

03 گس 

8 00 9005 0 7 0005 
(۷6) طاواد اد ۱۷۵۲۵ هت 
ج‌ ج‌ 
شکل ۸ مقايسة نتایج آزمایشگاهی کرنش مسطح دینامیکی شکل ٩‏ مقايسة نتایج آزمایشگاهی سه‌محوری دینامیکی ماس 
زهکشی‌شده ماسة متراکم با نتایج مدل ۳2 و مدل ۳2 اصلاح شده» متراکم با مدل ۳2 و مدل ۳2 اصلاح‌شده در شرایط زهکشی‌شده 
5-6 256» کرنش جانبی- کرنش قائم الف) آزمایش ماسودا و 2-5 0۵50 الف) آزمایش سه محوری دینامیکی ماسودا و همکاران 
قمکازان 22 تب هس سای با من 7 اصلاح‌شده و ج) یات 221 ب) مدل ۴2 اصلاح‌شده وج) مدل ۳2 
سازی با مدل ,۳۲7 
03 
2۵0-0 0 7701/6۷۲2 
جر عنلهی امتحمتر1 08 تس اب 
۷ 02 
2 ۳ 
3 


8 


هه 0125۱ > ی 


۶ ,51۲51۳ ۳۱۵۲۱2۵۳۵۵1 
ِ ه‌ِ 
4 ۳ 


ه‌ِ 


(۵0۵,۱ ۷۵ 1201۵)/( ۲-1۱ نی (قطط)م٩‏ 


۴8۵۶۵2059 
02 
06 04 02 00 02 0,4 
(90) > ,90۲210 اهه‌نا۲ع۱۷ 
0005 0 5 ۲ 
۱ (الف) 
(الف) 93 
88 ری 
۳۵0۵ ۳۸ ۱۱۵۵1۱۲۵0 و هروه ۳1 
23 
21 02 
مّ 
3 > 
كِ 
1.- 8 
3 ِ 
2 8 
1 3 
م 
1 ک 
ِ 5 
عِ ۰ 
۳ 
2 
ِ 0 
ب 
ت 
۳2 
4 
0005 0 005 ۹ 
, ,9 0۳ 04 02 0 02 04- 
1 1 2۳ 1 ۳ ۳ 
0 5۸۴ (۷0) دناد اجع۵ ۱۷۵۲ 
(ب) (ب) 
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۸ 
03 
02 
53 
۳ 
9 
تا 
ها 
8 
4 
۳9 
ت_ 
5 
۳ 
01 
05 04 02 0 02 رم 02 
(۵) 5۳2۱۳ ۱۷۵۲۱۵ 
‌ 


شکل ۱۰ مقایسُ نتایج آزمایشگاهی کرنش مسطح دینامیکی 
زهکشی‌شده ماس متراکم با نتایج مدل ۳2 و مدل ۳2 اصلاح‌شده 

2-5 806 کرنش جانبی - کرنش قائم. الف) آزمایش ماسودا و 
همکاران [22]. ب) شبیه‌سازی با مدل 2 اصلاح‌شده و ج) شبیه - 


سازی با مدل ,۲7 


اصلاح پیشنهادی دوم: اصلاح اثر رچتینگ 
رچتینگ به حالتی گفته می‌شود که در آن. نمونه تحت 
بارگذاری متناوب کرنش‌های خمیری پیش‌رونده ایجاد 
شود [24]. در مراحل اولیٌ وقوع زلزله معمولاً شتاب 
ناشی از زلزله پایین است و درنتیجه دامن کرنش در 
وهای تز تکیکی ‏ بانین ش ناگ فلع ی اب 
شرایط مقادیر کرنش خمیری تجمعی بالایی را نشان 
می‌دهد. به‌طوری که منجر به بروز تغییرشکل‌های قابل- 
توجه خواهد شد. تغییرشکل پیش‌بینی‌شده در این شرایط 
تشر از تغییرشکل وافعی است: شببه‌سازی. آزمایش‌هاق 
دینامیکی المانی برروی خاک‌های دانه‌ای توسط مدل ۳2 
این موضوع را نشان می‌دهد. 

پارامتر م13 برای اعمال اصلاح اثر پدیدة رچتینگ در 
مدل ۳2 معرفی می‌شود. این پارامتر به پارامترهای ,ولآ و 
بل به‌ترتیب در مرحلهٌ بارگذاری مجدد و باربرداری 
ضرب می‌شود. مقدار :۲۷ طوری تعریف می‌شود که تابع 
کرنش انحرافی خمیری تجمعی (]0ع4] [ < 5) باشد: 


(۲۲( 0 + (6و8-)۵۳ 02 < م11 


نشريه مهندسی عمران فردوسی 


معرفی دو اصلاح موثر برای بهبود پیش‌بینی رفتار دینامیکی حاک‌های... 


در رابطةٌ بالا 82 و 8 پارامترهای جدید مدل 
رفتاری هستند. ضریب خمیری بارگذاری مجدد و 


باربرداری در مدل رفتاری ۳2 به‌صورت زیر بازتعریف 


(۲۳( م«وا و + لاملا ع ۲ 
و 
۳ 
13 9 سا ْ () که 
جع 
(۲۵) 1 > ۳3 101 و0 ۳ < ۲ 


در سیکل‌های اوليةٌ بارگذاری و تغییرات محدود 
تلف ممقلاز گنت تخمیر ابخرآقی تیم یا زین 
است. بنابراین لو در نتیجه مقادیر و ن] بالا خواهد 
بود و تغییرشکل‌های برگشت‌ناپذیر نیز کوچک خواهند 
درنتیجه مقدار 11 کاهش پیدا می‌کند و نزدیک عدد یک 
می‌شود؛ بنابراین اثر پارامتر لا بر ضرایب خمیری 
بارگذاری و باربرداری کم می‌شود و درنهایت بسیار ناچیز 
می‌شود. پارامترهای 82 و وق طوری انتخاب می‌شوند که 
شیب نمودار تنش- کرنش در بارگذاری‌های بعدی 
بیشترین تطابق را با آزمایش داشته باشد و هم‌چنین تعداد 
سیکل‌های بارگذاری لازم در شبیه‌سازی و آزمایش برای 


تولید یک مقدار مشخص کرنش خمیری یکسان باشد. 


اعتبارسنجی اصلاح پيشنهادشده برای پدیده رجتینگ 
از یک آزمایش کرنش مسطح دینامیکی [25] و یک 
آزمایش سه‌محوری دینامیکی [26] برای صحت‌سنجی 
اصلاح پیشنهادی استفاده شده‌است. در ادامه آزمایش‌ها 
توسط هر دو مدل رفتاری ,۳2 و ۳2 اصلاح‌شده شبیه- 


سازی شده و با نتایج آزمایشگاهی مقایسه شده‌اند. 


شبیه‌سازی ازمایش کریمی و همکاران 
در اين بخش از یک آزمايش کرنش مسطح دینامیکی که 


سال سی و دوم شمارهة یک» ۱۳۹۸ 


امین ایرجی 


تویورا و در شرایط زهکشی‌شده انجام شده برای 
صحت‌سنجی اصلاحيهةٌ پیشنهادی استفاده می‌شود. 
مشخصات ماسه شامل: صه 202و 0.961عمرم 
01 زرم . 022.635 ۰ 9076 ارتفاع نمونة 
آزمایش 0 ۱ عرض آن حصحط 1۰ و طول آن ۸۰ 
می‌باشد. تدش جانبی ثابت اعمالی به نمونه 1۳2 14 است. 
آزمایش شمارة ۳۹002 (کریمی و همکاران [25]) برای 
شبیه‌سازی انتخاب شده است. 

جدول () پارامترهای مدل ۳2 را برای این آزمایش 
تضان می‌دهلند خیرات کرنتین: حجفی: هن ترایر. کرنسن 
محوری برای این آزمایش گزارش نشده‌است؛ بنابراین 
مقدار معا از سعی و خطا به‌دست آمده‌است. 

شکل (۱۱) نتایج آزمایش کرنش مسطح دینامیکی 
را نشان می‌دهد. همان‌طور که مشاهده می‌شود» پس از 
۰ سیکل بارگذاری و باربرداری, بارگذاری به‌صورت 
یک‌سویه تا مرحلة گسیختگی ادامه می‌یابد. شکل (۱۱- 
نب تقیان هی ده که خی ملل: مطاسین از ٩‏ سیک 
بارگذاری و باربرداری. گسیختگی روی می‌دهد. 
درحالی که بااستفاده از مدل ,۳2 اصلاح‌شده. ۵۰ سیکل 


۸۷ 


بارگذاری و باربرداری منجر به کرنش ۰ می‌شود که 


به نتایج آزماشی فزدیکاتر امتا: 


حدولعء پارامترهای مدل ۳7 و 2 اصلاح‌شده برای فایبنه متراکم 


تویورا در شبیه‌سازی آزمایش کرنش مسطح دینامیکی 


پارامتر مدل 7 
دول تم بابه(۳9 6 ۱/۰۹۹ 
مدول برشی پایه(۳۵) ,6 :2-۱3" 
پارامتر ثابت مدلم1۸ ۰/3۹ 
معط لت مان بل ۸ 
ضریب خمیری اولیه (بارگذاری) م11 ۱۳/۹ 
پارامتر ثابت مدل (باربرداری) من[ وتم ال 
پارامتر ثابت مدل (بارگذاری مجدد) *۷ ۰/- 
پارامتر ثابت مدل (بارگذاری مجدد) ۷۲ ۱/۱ 
پارامتر ثابت مدل (باربرداری)ن ۷ #7« 
پارامتز خایت مدل م8 0/۰ 
پاژاشد اقایت‌تمدلرق ۲/: 
پارامتر ثابت مدل (بارگذاری و باربرداری) 8 ۳33/۰ 
پارامتر ثابت مدل (بارگذاری و باربرداری) 83 ۳/۰ 


* فقط مدل 7 *فقط مدل ۳7 اصلاح شده 


سس ولووا م۵۵۵۵ آن ۵۷ وصالهه و۷۵۵ ۱۰ آن 669۱00 


(:6۷۵۱6 50) رز + 


880۷01 )20 6۷۵۱۵6( 


(90002 


"۰:01 


(۹60002 1۵۲) 24 رس 


96( ]6*۷16۲0۵۵۱ ۱۲۵۱90۵66۲[ 


(الف) 
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شکل ۱ آزمایش کرنش مسطح دینامیکی زهکشی‌شده برروی ماس تویوراه نمودار تنش انحرافی- کرنش محوری: (الف) آزمایش کریمی و 
همکاران 1251 (ب) مدل 72و (ج) مدل 7 اصلاح شده 


شبیه‌سازی آزمایش لینگ و همکاران 

لینگ و همکاران [26] یک آزمایش سه‌محوری دینامیکی 
زهکشی‌شله برروی ماس توکاچی (۲0100) انجام 
دادند. در این آزمایش تنش جانبی 1۳2 ۷۰ و دانسيتة 
ذسبی نمونه آزمايش ۵۵ است. جدول (۵) پارامترهای 
مدل ,۳2 و ۳2 اصلاح شده را برای ماسة توکاچی ذشان 
می‌دهد. این پارامترها براساس نتایج آزمایش سه‌محوری 
تزورون ماه کالت هلان 

شکل (۱۲) رفتار ماس توکاچی را در آزمایش سه- 
محوری نشان می‌دهد. نمودارها شامل تغییرات تنش 
انحرافی دربرابر کرنش محوری و کرنش حجمی دربرابر 
کرنش محوری هستند. شکل (۱۳) شبیه‌سازی آزمایش 
سه‌محوری با مدل 2 را نشان می‌دهد. هم‌چنین شکل 
(۱۶) نتایج شبیه‌سازی آزمایش سه‌محوری را با مدل ۳7 
اصلاح‌شده نشان می‌دهد. در شکل (۱۲-لف) مقدار 
کرنش محوری پس از ۵ سیکل بارگذاری و باربرداری 
یا 

مدل ,۳2 1 کرنش محوری و مدل ,۳2 اصلاح‌شده 
کرنش محوری پس از ۵ سیکل بارگذاری و 
باربرداری را محاسبه می‌کند. کرنش محوری پیش- 
بینی‌شده توسط مدل 32 اصلاح‌شده به‌میزان واقعی آن در 


نشريه مهندسی عمران فردوسی 


آزمایش نزدیک‌تر است. کرنش حجمی تجمعی پیش- 
بینی‌شده با مدل ,۴2 اصلاح‌شده در شکل (۱۶-ب) برابر 
است که تطابق مناسبی با مقدار واقعی آن در 
آزمايش یعنی /۱/۷ دارد. کرنش حجمی تجمعی پیش- 
بینی‌شده با مدل ,2 معادل /۳/۵ است. 


حدول ۵ پارامترهای مدل 2 و مدل ۳7 اصلاح شده برای شبیه - 


سازی آزمایش سه‌محوری دینامیکی برروی ماس توکاچی 


بار اف ما مقدار 
مدول حجمی پایه(0۳ وک 2 
مدول برشی پایه(۳۵) » وت 
پارایع ایک سول 9۳ 
شیب خحط حالت بحرانی و ۱/۳۵ 
ضریب خمیری اولیه (بارگذاری) م1 ۹ 
پارامتر ثابت مدل (باربرداری) من 11 ۱/۳ 
پارامتر ثابت مدل (بارگذاری مجدد) 7۹ ۱1/۰ 
پارامتر ثابت مدل (بارگذاری مجدد) 7 ۳/۰ 
پارامتر ثابت مدل (باربرداری)ی 7۷ ۷/۰ 
بازامتر کات مد وق ۱/۰ 
تار قایت مان ۶ ۰/۱ 
پارامتر ثابت مدل «بارگذاری و باربرداری) 84 ره 
پارامتر ثابت مدل (بارگذاری و باربرداری) 64 ۵/ 


* فقط مدل ۳2 "فقط مدل ۳2 اصلاح‌شده 
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زهکشی‌شدهة لینگ و همکاران [26]: (الف) تغییرات تنش انحرافی (الف) 


در برابر کرنش محوری و (ب) تغییرات کرنش حجمی دربرابر 


کرنش محوری 
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)- 

ت 
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5 

(ب) 
و شکل ۱۶ شبیه‌سازی آزمایش سه‌محوری دینامیکی لینگ و 
ا ا را ا را 0 
«نووه ازجم همکاران [26] با مدل ۳2 اصلاح‌شده: (الف) تغییرات تنش انحرافی 


در برابر کرنش محوری و (ب) تغییرات کرنش حجمی دربرابر 
(الف) 7 
کرنش محوری 
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نتیجه گیری 

دو اصلاح مزثر برای بهبود پیش‌بینی رفتار دینامیکی 
خاک‌های دانه‌ای متراکم و هم‌چنین بهبود پیش‌بینی 
کرنش‌های خمیری پیش‌رونده تحت بارهای متناوب 
(رچتینگ) ارائه شد. در ابتدا دو سری آزمایش سه‌محوری 
دینامیکی برروی ماسهٌ سست و متراکم انتخاب شد و 
رفتار آنها توسط مدل رفتاری ۳2 پیش‌بینی گردید. در 
طول این بررسی مشخص شد که مدل رفتاری ۳2 قادر 
به شبیه‌سازی پدیده سست‌شوندگی ماسة متراکم در طول 
بارگذاری متناوب نیست و لازم است یک اصلاح موثر 
بدین منظور در مدل اعمال شود. اصلاح مذ کور با شبیه- 
سازی دو سری آزمایش سه‌محوری دینامیکی برروی 
ماسه صحت‌سنجی شد. 

مدل ۳2 در خاک‌های دانه‌ای تحت بارهای متناوب 
و با سیکل‌های تنش نسبتاً کوچک کرنش‌های برشی 
حمیری پیش‌رونده بیشتری را نسبت به آزمایش پیش‌بینی 
می‌کند. این مسئله ناشی از نقص مدل ۳2 در پیش‌بینی 
مطلوب پدیدةٌ رچتینگ به خصوص در سیکل‌های 


2 


معرفی دو اصلاح موثر برای بهبود پیش‌بینی رفتار دینامیکی حاک‌های... 


کوچک بارگذاری و باربرداری است. برای بهبود پیش - 


اصلاح مذ کور نیاز به دو پارامتر جدید دارد که منجر 
به بهبود پیش‌بینی پدیدهٌ رچتینگ می‌شود. ابتدا از یک 
آومانقن سه‌محوری دینامیکی فرضی استفاده شد. سپس 
از یک آزمایش کرنش مسطح دینامیکی و یک آزمایش 
سه‌محوری دینامیکی برای این صحت اصلاح اعمال‌شده 
استفاده سل در خصوص اصلاح دوم بایستی اشاره نمود 
که اصلاح مذکور قادر به پیش‌بینی دقیق پدیدة 
هیسترزیس در سیکل‌های بارگذاری و باربرداری 

از مدل رفتاری اصلاح شده می‌توان برای شبه‌سازی 
دینامیکی سازه‌هایی استفاده کرد که عمدتاً از مصالح دانه- 
ای متراکم تشکیل شده‌اند. هم‌چنین در سازه‌هایی که 
تحت بارهای دینامیکی نسبتاً طولانی قرار گرفته و 
کرنش‌های خمیری جمح‌شونده در آنها قابل توجه است؛ 
می‌توان از مدل اصلاح شده استفاده کرد. 
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(2001) ,336 

۵ لوصح عصتاام۷۲۵0 ممتنحعللهما صتقتاک ریم راهم اصم2 عک ال رتمافه۲ وبط معصفط2 ر.ظ اههد 
(1998) ,266-280 ۲0.۰ 3 0۰ ,22 ,۷۵۱ روعنصهط)ع]۱۷ آمممتاه) تامهم و وعام‌هه 4صقه ٩200۳2660‏ صد معتاوومزظ 
عمط مصحو ۶ ممرتحطفظ ممتاهعز/نفوهز مج پمصمل‌مممعز آمع1-عناووعزظ۳. وبا بتانا 200 .1 ,عط] 
,851-860 .30 و8 .۱0 ,129 ۱۷۵۱۰ ۱۱۵6۱۵ و۵8۵۵ زن 0۵۵ ,۱۷۲۵۵۵ اته‌تاد۱۵ظ۳ 0۵۵۵1۵11760 
2003(۰) 

۳۱2601010 0۵۵۵۵11760 )۳1 ممممهنا مصتااه0م]۱۷ ر.ظ بهت۷ ,۷ رتماعهط ریگ متطتصعقحه ] ببظ و.ل ,1۷۲۵۲000 
(2005) ,۴0.32-42 ,4 ,۷۵۱ بعتصععاهع06 تنل مصمنلمآ ماوزنتن؟۴] 

۲ صطاتمواهه امئمم۱۲0 060۵21760 ۲۸۵ ربظ ,۱۷۵۲۵00 2عصمصع۲ صه رب رتمافهط ۷ بتصجم] ر.ظ بهتذآ۷( 
۰ ,77 ۷۵۱۰ ,۱۱۵8 ۲۵۱ ۷۵۱۵۲ ل ۱ 52905 ۶۵۲ ۱۷۵0۵1 انامه[ 062۵2112760 ۵ 0۴ طمتهتععاصا ما 
(2009) ,1201-1223 ۳0۰ ,9 

۴ ۵)2000 مجمحجامتامی عمط ممتاعع ]عاونا مصجو رم م۱۷۵0 7۵تاتاتاعجمن ۸ بب۷ رتکتاصضا ک بط راتهاطوم ] 
2009(۰) ,215-223 ,2 ۱0۰ ,36 ۷۵ م۵0 هک ۱۷۵۵۱۵۵ ,ومع وومتوو آهمامطز:۳۲ 
متاعجام»توهاط؟ طز م۷۲ مصع عاعله لعفمتاتیت) ۸ به. بتتقلطاعقا ک بض مامت 
2014(۰) ,2807-2825 ,15-16 ۱۱0۰ و51 ,۷۵۱ ۳0اه آصه ولآ هن مرول ۱۱۹۵۵۵۵ عمجم 
ممنهاصصته متصصجصو ۱۱ عصهو مقصونا م] صم۱۷]۵01)162 ۸ ر.ظ رهتصتصمووم۲۱ 907601 220 .۵ رطمصه۳2۲2 رنه رازه1[ 
2014(۰) ,343-353 .۲۵ 2 .۱0 ,9 ۷۵۱ بهعنصطهعا۵ع0 مامظ ,۱۷۱۵061 مقط)-کم۲۵۵)۵۲-10۵110۱ ۶و ندومن 
۷۱0۲۵6 ۷۵ ۵ع ص52 0۶ ۱۷۱۵061 ع7تانتاعدمی وطتم‌تاعها۲ 060۵۳211760 ۰۸ رب رل‌تصه۲۱ مه رب رتصحتهم0) 
.(2015) ,133-145 .0 و2 .۱0 ,13 ۷۵۱۰ ,1۱۵۱۱668 ۷ ۵۴ ول ۲۱۱12۵۵ 

۲ ۱۷۱۵06 ۷۵)نتاوصم) المتافهاظ 2۵0ناهتمهع0 وه ۵۶ امممصمم۰۱۵۷۵۱ ره رتتصه۲۱ 24 ر.ط مطفطلهماصه۲۷ ۲۴۵2 
و 6060۱۵۵۱۱۵۵ ]0۵ 0۵ ۳۱۵۳۱۱۵۵۵ 0۵۵۵۵۵۲ مها لهعتانت عطلونا فلصهو 0ماصمصمن 
2013(۰) ,364-373 ۳۵۰ ,4 .0 ,۲۷۵۱۰7 

2 ,۸۵۵۲0200۲ ۱۷۲۵00 ع۷تابمتاعجم 2 مصمتامد‌]مدونا ۳۱۵ عصتام‌نلع:۳. ,۸ رتلققا مه .ده مقصوزعل 
2017(۰) ,793-808 ۴۵۰ ,4 .0 ,12 ۷۵۱۰ بعتصطه‌ماهع0) 

۰ 0۲۵۲۵۹5101) أمتیمز: ۲ مد 0صحه نم م۷۲۵0 عتاناتاوممن هر رب .1۲ ,۷۷۵۵۵ 20 ریک ب0۷2 


.(1979) ,295-278 ۲0۰ و3 .۱۱0 و3 ۷۵۱۰ 60۵60 ۱ و۱۱۵۵ م0۵ وه 0۱۵۵۲۱۵۱۵۵۵ 
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0 924 ۵۲ عممتاقای؟ «مصداهانل توعد آماممصتنمم ۳۶ .۷۲ رماقه و.۲ بق‌امتاعه] پگ .ظ .۲ مصقطل۳:۵2 
.(1989) ,45-64 ,۴0 م1 .0 ,29 ,۷۵۱ رعطمتا02ع۳۵0 لمح واتمو ر عصتلعما متلهن ما 

صنععگ مصهاط طاً فصق ۵۶ عمتقطهظ صتقت ای وب رماع عک ویک م2۳30202 ۷ .۳ م2050101] ب.] ,02امع]۷[ 
,31-45 .۲0 ,5 ,0 و39 ,۷۵۱ ,کعمتعهمت۲۵ مج ولزمو , عاوع1 عمتجم متام 4صه ممتقمم ۲ بممتوی‌تموصمن 
.(1999) 

م۵ موم متل0۵ فصن فصعی ۵۶ ممتاممصهم]ع(1 0مصنعنا. _ ویک رهتقطتطع فصه و۲ رقامتاواع ]1 
19740(۰) ,51-65 .۲0 ,3 .0 ,14 ,۷۵۱ رقصمناهصت۲0 مه ولژو ر "صمتامع:ز(] 

۵0826 صنجگ متافهاظ! از عاتمو فوماممتوعطمت 1۵۲ ۷۵061 متافهام۰۳۱۷۵۵ ربا رقاتع۳۱ 20 بیظ رقتصتاصهمزل 
٩۲0. 4, ۳0۰ 279-299, )1997(,‏ ,2 ۷۵۰ ,علهتتم)۱2 لمصمتامن -ع۷تومطمن ۵ ومنصهطمع]۷( 

۵ 8 ۸217516 معفقصصا لاه فاوع1 ممتوومتمموم صتجتای مصعاط. ر.ژ رتعله‌عمک عک ر.] ,9260 ول متطاته 
2005(۰) ,81-92 .0ظ ,38 .۲0 ,هل ]0 ۱0/۱۵1 9200 ما100 

حالنا۲ ۵۶ عصمتحانهته صمجصماهمانص۳ رت رفظ 4ص رب بتاشا رثا ریافمتط‌طوما ری رعط۷2 وی .1 رع19] 
۵9۱۵ ,عصنلدصا ملمتا مج عم ولله ۷۷ ممتمتجای؟ ارم 0عمعمصنم؟ ‏ 1مم۱ه-تقلن0م] ۷( قلوعو 
2010(۰) ,653-661 .۴0 ,5 .۱0 ,136 ۷۵1 ,۱۸۵۵۱۱۵ ۲۱۵۱۱۱6۵۲1۱8 
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